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We have mapped from the quantum to the classical limit the spin excitation spectrum of the 
antiferromagnetic spin-1 Heisenberg chain system CSN1CI3 in its paramagnetic phase from T = 5 
to 200 K. Neutron scattering shows that the excitations are resonant and dispersive up to at least 
T = 70 K, but broaden considerably with increasing temperature. The dispersion flattens out with 
increasing temperature as the resonance energy A at the antiferromagnetic wave-vector increases 
and the maximum in the dispersion decreases. The correlation length £ between T = 12 and 50 K 
is in agreement with quantum Monte Carlo calculations. £ is also consistent with the single mode 
approximation, suggesting that the excitations are short-lived single particle excitations. Below 
T — 12 K where three-dimensional spin correlations are important, £ is shorter than predicted and 
the experiment is not consistent with the random phase approximation for coupled quantum chains. 
At T = 200 K, the structure factor and second energy moment of the excitation spectrum are in 
excellent agreement with the high-temperature series expansion. 



PACS numbers: 75.25.+Z, 75.10.Jm, 75.40. Gb 
I. INTRODUCTION 

Low-dimensional antiferromagnets are ideal systems in 
which to investigate the effect of strong quantum fluctua- 
tions on the dynamics of interacting spin systems. Quan- 
tum fluctuations are similar to temperature fluctuations 
in three-dimensional (3D) magnets in that they disor- 
der the spin system and prevent magnetic order. Like 
the paramagnetic phase of a 3D magnet, the quantum- 
disordered phase of a low-dimensional antiferromagnet 
has unbroken translational symmetry, but while the ex- 
citation spectrum of a 3D paramagnet consists of broad 
magnetic excitations, the elementary magnetic excita- 
tions of a low-dimensional aajtiferromagnet are spin-1/2 
or long-lived spin-1 particles lit] 

One-dimensional (ID) antiferromagnets have been 
studied for at least the last 30 years. They were first per- 
ceived as simpler model systems than their 3D counter- 
parts, but it is now clear that their low-temperature spin 
dynamics is far more complex. The quantum-disordered 
ground state of ID antiferromagnets is usually a highly 
correlated spin singlet. Haldane, in his seminal 1983 pa- 
per, conjectured that the spin dynamics of ID antiferrcp. 
magnets crucially depends on the spin quantum number :□ 
The half-integer spin Heisenberg antiferromagnet has a 
gapless excitation spectrum while integer-spin chains ex- 
hibit an energy gap. Haldane's conjecture of a spin-1 
gap was then, confirmed by numericalB'tl and experimen- 
tal studiesacl and is now generally accepted. 

The low-temperature quantum-disordered ground 



state of the antiferromagnetic spin-1/2 chain is "almost 
ordered" and its correlations fall off with distance as a 
power law. The elementary excitations are spinons car- 
rying spin-1/2 which correspond to moving domain walls 
between regions of coherent quantum states. In contrast, 
the correlations in antiferromagnetic spin-1 chains fall off 
exponentially with a correlation length of about 6 sites at 
zero temperature. The elementary excitations are mas- 
sive triply-degenerate spin-1 particles, called Haldane ex- 
citations, which can be pictured as domain walls break- 
ing theJiidden spin string order that is inherent in spin-1 
chains.B The triplet excitations follow the relativiatic dis- 
persion of the non-linear sigma model (NLctM).Cj 

The temperature dependence of the excitation spec- 
trum of antiferromagnetic chains has been measured 
mostly for T < J, where J is the exchange interac- 
tion between the spins along the chain direction. For 
KCUF3, an antiferromagnetic spin-1/2 chain system with 
J = 17 meV, inelastic neutron scattering up to T = 
200 K ~ J showed that the low-energy excitations are 
broadened at high temperatures bui_can still be described 
by a low-temperature field theory]!!] 

For spin-1 chains, the temperature dependepcje,|-af 
S(Q,uA has beep-, investigated fot CsNiCyMliy, 
NENPB, AgVPaSgN and YaBaNiOgtil for temperatures 
T up to at most T ~ J . The temperature dependence 
of the excitation for T < 0.3 J was studied for AgVP2S6 
and Y 2 BaNi0 5 (J = 60 meV and J = 24.1 meV, respec- 
tively). These measurements give the behavior of the 
low-temperature quantum-disordered phase. A further 
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inelastic neutron scattering study of AgVP2S6 showed 
that the staggered susceptibility decreases, with increas- 
ing temperature up to T = 200 K ~ 0.3 JS3 A single crys- 
tal experiment on ISBaNiO.-, showed that the energy gap 
observed by neutron scattering is renormalized upward 
with increasing temperature up to T = 80 K ~ 0.3 J.E3 

CSN1CI3 and NENP have much smaller exchange in- 
teractions, J = 2.28 meV and J — A meV, respectively, 
and were studied over a wider temperature range. In 
NENP the upward renormalization of the energy gap i-u, 
the spectrum was measured up to T — 20 K ~ 0.43J,EZl 
and the temperature dependence of the correlation-length 
has been determined up to T = 100 K ~ 2JJB The 
temperature dependence of the gap in CsNiCbj at the 
3D ordering wave-ijeator was previously studied up to 
. 75 ,7" vrrH and at the ID wave- vector up to 
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Surprisingly a detailed study of the excitation spec- 
trum of antiferromagnetic spin chains for T > J and 
in the high-temperature limit has not been made. We 
showed recently that the Haldane gap in CsNiCl3 per- 
sists m a resonant feature up to at least T = 70 K or 
2.7J.E3 Here we present more extensive measurements 
at low and high temperatures including a comprehensive 
study of the spin dynamics of the ID magnet for tempera- 
tures T > J. CSN1CI3 has an exchange interaction J that 
makes it an ideal system for studying high-temperature 
fluctuations, not only because the whole magnetic exci- 
tation spectrum can be measured with thermal neutrons 
but also because temperatures between J and ~ 10J 
(room temperature) can readily be accessed. 

Section II of the paper describes magnetic properties 
of CSN1CI3 and gives the configuration of the neutron 
scattering spectrometers. Section III gives a.-thorough 
analysis of the measurements presented earlier. E3 The in- 
elastic time-of-flight neutron scattering measurements of 
the excitation spectrum between T — 12 and 50 K are 
reported in Section IV. In this temperature range, reso- 
nant and dispersive excitations were observed for a wide 
range of wave- vector transfers along the chain. In Section 
V, we describe the measurement of the dynamic struc- 
ture factor, S(Q,uj), of the paramagnetic scattering at 
T = 200 K. The results are summarized, discussed and 
put into an overall context in Section VI. 



II. EXPERIMENTAL DETAILS 

CsNiChj is one of the best examples of a gapped 
weakly-coupled spin-1 chain system. The crystal 
structure is hexagonal, Df h space group, with low- 
temperature lattice constants a = 7.14 A and c = 
5.90 A. The exchange interaction along the c-axis, J = 
2.28 meV,Ej is much stronger.-t.han the interaction in the 
basal plane, J' = 0.044 meVj3'E3 making CsNiCl 3 a sys- 
tem of weakly-coupled spin-1 chains. The Hamiltonian 
can be written as 



The weak Ising anisotropy D = 4 fieV is small enough 
that CsNiCl3 is a good example of an isotropic Heisen- 
berg antiferromagnet. 

Because of the interchain couplings CsNiO.3 orders an- 
tiferromagnetically below T N = 4.84 K.E3 Above T N , 
CSN1CI3 displays the characteristics of an isolated anti- 
ferromagnetic spin-1 chain: the excitation spectrum ex- 
hibits the Haldane gap, the correlations fall off exponen- 
tially and the magnetic susceptibility shows a broad max- 
imum near 30 KH3 

The single crystal of CsNiCLj 5 mm x 5 mm x 
20 mm was sealed in an aluminum can containing he- 
lium gas to prevent absorption of water. The experiments 
were performed using the cold-neutron RITA triple axis 
spectrometeio at the DR3 reactor of the Ris0 National 
Laboratory, Denmark, and the chopper time-of-flight 
spectrometer MARI at the pulsed spallation source ISIS 
of the Rutherford Appleton Laboratory, United King- 
dom. 

For the experiment using the RITA spectrometer, 
the sample was mounted in a cryostat with its (hhl) 
crystallographic plane in the horizontal scattering plane 
of the neutron spectrometer, and the temperature was 
controlled to an accuracy of 0.1 K between 1.5 and 
70 K. The energy of the neutrons from a cold source 
was selected by a vertically focusing pyrolytic graphite 
monochromator. A rotating velocity selector before 
the monochromator suppressed unwanted neutrons that 
would be Bragg reflected by its higher order planes. Su- 
permirror guides with n = 3 (6 C = 1.2° at 4 A) were 
located before and after the monochromator, and the 
beam was 20 mm wide. The scattered beam was filtered 
through cooled beryllium, passed through a 50 mm wide 
1° Soller collimator, and was analyzed by reflection from 
the central two blades of a 7-component flat pyrolytic 
graphite analyzer aligned so that each blade reflected 
the same energy of neutrons, 5 meV. The Soller geome- 
try meant that the beam reflected by the analyzers was 
about 20 mm wide and located in the central 30 strips of a 
position sensitive detector. Because the detector was 120 
pixels wide (~ 120 mm) the side strips were used to es- 
timate and subtract the temperature independent back- 
ground. Turning of either analyzer or monochromator 
from their reflecting position confirmed that the side pix- 
els gave a good representation of the background arriving 
at the centre of the detector. With this arrangement the 
low-energy resolution was typically 0.26 meV (FWHM) 
in energy. The longitudinal and transverse wave-vector 
width of the (002)-peak was 0.01 and 0.016 (FWHM in 
reciprocal lattice units) and the calculated vertical reso- 
lution was 0.22 A -1 . 

For the experiment using the MARI spectrometer the 
sample was oriented with its (hhl) plane in the vertical 
scattering plane, and the c-axis was either parallel or per- 
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pendicular to the incident beam direction. The incident 
neutron energy was either 20 or 30 meV with the Fermi 
chopper spinning at 150 Hz. The energy resolution was 
then 0.35 or 0.4 meV, respectively, as determined from 
the full width at half maximum (FWHM) of the quasi- 
elastic peak. The resolution in wave-vector transfer at 
zero energy transfer was typically 0.02 A -1 along both 
the c* and the [110] directions and as large as 0.19A -1 
perpendicular to the scattering plane if only the central 
detector bank was used. Both the energy and the wave- 
vector resolution improved with increasing energy trans- 
fer. The measurements were performed with the sam- 
ple at temperatures between T = 6.2 and 200 K. The 
scattering was measured for a total proton charge be- 
tween 4000 and 8600 /iAh at an average proton current 
of 170 nA. 

MARI has three detector banks two of which are out 
of the vertical scattering plane. For T < 25 K only data 
recorded by the central detector bank was included in 
the analysis because the two side banks sample an out- 
of-plane wave- vector transfer where the neutron scatter- 
ing is different. All three detector banks were used for 
the measurements performed at temperatures T > 50 K 
where the 3D dispersion is small and good resolution per- 
pendicular to the scattering plane is unnecessary. 



trum at T = 6 K that an intrinsic Lorentzian line-shape 
gives a better description of the, measured data, as pre- 
dicted by Damle and Sachdev.E-3 
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FIG. 1. Scattering from the Haldane excitation at 
Q = (0.81, 0.81, 1) for T = 6 and 9 K as a function of energy 
transfer. The data were measured using the triple-axis spec- 
trometer RITA as described in the text. The intensity was fit- 
ted to a Lorentzian (left) and a Gaussian-Xright) cross-section 
convoluted with the resolution functionO 



III. EXPERIMENTAL RESULTS 
A. Excitation at Q c = 1 

Close to the antiferromagnetic point, denoted Q c = 1 
because the spin spacing is c/2 along the chain, the ex- 
citation spectrum of an antiferromagnetic spin-1 chain 
in its quantum-disordered phase is dominated by well- 
defined, gapped spin-1 excitations (Haldane excitations) 
which satisfy a relativistic dispersion relation. We 
showed in a recent publicationESl that for T < A, the 
Haldane excitation at Q c = 1 exhibited a Lorentzian 
spectral form. We showed that its life-time agrees well 
with a semi-classical theoryO for gapped systems pro- 
vided we replaced the classical dispersion by the more 
accurate relativistic dispersion of the NLerM. At higher 
temperatures the Haldane excitation evolves into a broad 
resonant feature. In this paper a more detailed analysis 
of the data is presented and results based on two dif- 
ferent spectral functions are presented for temperatures 
between T = 5 and T = 70 K = 2.7J. 

The Haldane excitation was measured at Q — 
(0.81,0.81,1) where the Fourier transform of the inter- 
chain coupling vanishes, and within a random phase ap- 
proximation (RPA) the chains behave as if they were 
uncoupled.H In Fig. [j], the results are shown at T = 6 
and 9 K. The spectra were fitted to both an antisym- 
metrized Lorentzian and an antisymmetrized Gaussian 
line-shape, weighted by the Pose factor, and convoluted 
with the resolution function.c3 It is clear from the spec- 



The Lorentzian fits were performed with a double 
Lorentzian form, namely: 



S(Q,u) = A ■ (n(u) + 1) 

r r 



e(Q)) 2 +r 2 ( w + e (Q))2+r 2 



(2) 



where u) is the energy transfer and the Bose factor is 
defined by 



nH + l = (l-cxp(--^))- 1 



(3) 



The energy width is given by Y (half Lorentzian width at 
half maximum) and the energy of the excitations by e(Q). 
A is an overall scaling factor which scales as l/e(Q) for 
antiferromagnetic spin waves and also for the excitations 
of the NLtrM. With increasing temperature, the excita- 
tion broadens as may be seen at T = 9 K, where the 
line-shape can be described equally well by a Gaussian 
or Lorentzian line-shape. 

Another spectral form that describes damped excita- 
tions is the damped harmonic oscillator (DHO) function 
which can be written as: 



S(Q,u) 



A- (n(uj) + l)4Ye(Q)uj 

^ 2 -£DHo(Q) 2 ) 2 + 4r2c 



(4) 



£dho(Q) is the DHO energy and the other symbols are 
as in Eq. |^. The DHO function is the same function as 
the antisymmetrized Lorentzian defined in Eq. if 
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e DHO (Q) 2 = e(Q) 2 + r 2 



(5) 



For the same spectral distribution, the DHO energy 
£dho(Q) is therefore larger than the Lorentzian energy 
e(Q). 
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FIG. 2. Neutron scattering intensity at Q = (0.81,0.81, 1) 
at T — 25 and 50 K measured using RIFA. Fhe data is shown 
as a function of neutron energy transfer. Fhe solid lines cor- 
respond to a Lorentzian cross-section convoluted with the res- 
olution ellipsoid and fitted to the data. 



In Fig. y we show that neutron scattering at Q = 
(0.81, 0.81, 1) at T = 25 and 50 K is well described by a 
Lorentzian cross-section. The background was assumed 
to be constant and a Gaussian was used to account for the 
incoherent elastic scattering. The cross-section (Eq. ||), 
folded with the resolution function,E£l was fitted to the 
data to derive e(Q), T and A. 

The excitation energy e(Q) and width T are shown in 
Figs. H and [| as a function of temperature T. The Hal- 
dane gap energy, which at T = 5 K is 1.2 meV, increases 
rapidly between T = 12 and 30 K as the excitation de- 
velops into a broad resonant feature. At about 30 K, 
a temperature of the order of the exchange coupling J 
along the chain, there is a cross-over to a qualitatively 
different behavior and for T > 30 K the energy increases 
only slowly with increasing temperature. 

The observed gap energies are comparedJn Fig. ^ (solid 
line) with the theory of Jolicceur et alE3 based on the 
NLctM, calculated with the T = K gap energy of Ao = 
10.6 K we obtain from the exchange constant as Ao = 
0.406 J. The gap equation of NLctM is a statement that 
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is a constant in temperature where 



e(q) = y/A(T)* + vW- 



(6) 



(7) 



and use to obtain g at T = K. For T < A = 11 K 
(see inset), where the theory is expected to be valid, the 
observed energy matches the theoretical prediction ex- 
cept below 7 K where the gap energy is higher. The 
good agreement for 7 K < T < 11 K shows that the 
NLctM gives a good description of the low-temperature 
quantum-disordered phase for T < A. The disagreement 
below T < 7 K is likely due to 3D interactions,EII These 
give rise to low-energy modes as Tjv is approached whose 
thermal population will drive up the gap energy in the 
NLctM equation (Eq. g). 
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This implicit equation for A(T) can be used to calculate 
the gap relative to Aq which we take from experiment 
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FIG. 3. Excitation energy e(Q) for Q = (0.81, 0.81, 1)) as a 
function of temperature between T = 5 and 70 K. The ener- 
gies were obtained from fits to a resolution-folded Lorentzian 
or damped harmonic oscillator function. The solid line shows 
the upward ^normalization of the gap energy obtained from 
the NLerM.Ej The inset shows the low-temperature range 
where NLctM is valid. 



The measured widths T are shown in Fig. |j. The width 
is always smaller than the excitation energy up to at least 
T = 70 K because the upward renormalization of the ex- 
citation energy is fast enough tp-pvercome the increasing 
damping. Damle and Sachdevo have developed a the- 
ory of the damping of excitations in gapped spin systems 
which is valid as long T < A. In CsNiCl3 this tempera- 
ture range is limited because A = 0.41J = 11K and 3D 
effects are important close to the ordering temperature 
Tjv- As shown in Fig. ^ their theory, which approximates 
the NLctM dispersion by a classical form, A + cq 2 , pre- 
dicts considerably longer excitation life-times than those 
observed experimentally. However, we find that Damle 
and Sachdev's theory agrees well with the measured ex- 
citation width between T = 9 and 12 K if we replace the 
classical dispersion of the elementary excitation particles 
by the more accurate relativistic dispersion of the NLctM 
(Eq. |). Details of our calculations are given in Section 
MA. Below T = 8K, the observed width is consistently 
higher than that predicted by the theory. Probably this is 
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due to the 3D spin correlations which are most important 
close to Tjv- We expect that the exponentially greater 
thermal population of the low-lying 3D gap modes will 
shorten the collision life time relative to that for the gap 
of an isolated chain. 
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FIG. 4. Half width F, corrected for resolution, for 
Q = (0.81,0.81,1) between T = 5 and 70 K. The widths 
were obtained from fits to a Lorentzian or a DHO function 
and averaging the widths. The solid line is the prediction of 
Damle and Sachdev. The dashed line is our modification of 
their model to include the relativistic dispersion of the NLctM. 
The inset shows the low-temperature range up to T = 15 K 
where the NLcrM is valid. 



Above about T = 20 K ~ J the temperature depen- 
dence of the excitation width decreases, similar to that 
of the energy. This behavior probably arises because the 
lattice momentum cut-off is not included in. the theoret- 
ical model. The cut-off is infinite in Ref.cfl, and in our 
relativistic calucation it was limit to q = ir (see Section 
VA). An extension of the theory to higher temperatures 
would be valuable. 

Fig. H shows the integrated intensity of the Haldane 
excitation, S H (Q), for Q = (0.81,0.81, 1), the ID point, 
as a function of temperature (solid circles). The Sh was 
determined from the observed neutron scattering inten- 
sity by integrating the area under the Lorentzian cross- 
section for positive and negative energies. It includes all 
the scattering except below ~ 12 K where the spectrum 
contains a high energy continuum of 12(2)% in weight .c3 
Sh increases with decreasing temperature as the antifer- 
romagnetic fluctuations become stronger and the spec- 
tral weight is displaced to lower energy. The increase 
between T = 50 and 20 K is associated with the growth 
of the ID correlation length. However, at 10 K the tem- 
perature dependence of Sh has a peak and then decreases 
with decreasing temperature. This is correlated with the 
hold-up in the gap decrease noted for Fig. |L 

The structure factor S(Q) for Q = (0.81,0.81,1), 



which includes the continuum, has a lesser decrease at 
low temperatures (Fig. ||). S(Q) was determined by nu- 
merical integration of the observed neutron spectrum at 
positive energy transfers and adding intensity at nega- 
tive energy transfer calculated using the detailed balance 
relation. S(Q) and Sn(Q) are identical at high tempera- 
tures, but for T < 10 K , where the continuum scattering 
was observed, S(Q) is larger than Sn(Q) and approxi- 
mately constant. The increase in S(Q) on cooling is suit 
stantially less than that calculated for isolated chains .E£l 
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FIG. 5. Structure factor at the ID point Q = (0.81, 0.81, 1) 
as a function of temperatures. The solid circles show the 
integrated intensity of the Haldane excitation 5h(Q). The 
open circles represent S(Q) determined from the entire neu- 
tron scattering intensity. Inset: The susceptibility of the 
Haldane excitations for positive energy transfer, defined as 
Xh(Q) = J °° du}S K {Q,u>)/{n{e{Q)) + 1), multiplied with the 
excitation energy e(Q) shows a transition to a qualitatively 
different behavior at ~ 15 K. 



To test whether the spectral weight scales with \/e(Q) 
we divided Sh by the population factor n(e(Q)) + 1 
to obtain the susceptibility for positive energy transfer 
Xh and multiplied it with the excitation energy e(Q). 
This approximately removes the thermal factor to a de- 
gree that improves as the temperature and width de- 
crease. Our scaled function (Fig. ||)is a constant for 
T > 15 K, showing that the intensity scales with l/e(Q). 
For T < 15 K, however, xh is lower than expected from 
a l/e(Q)-dependence of the intensity, thus the spectral 
weight grows more slowly than l/e(Q). 

The reason for the behavior below-J 1 = 10 K is at 
present unclear. In Kenzelmann et al£3 we have argued 
that the continuum scattering is generated by the in- 
terchain coupling which has a strong effect on the dis- 
persion for temperatures close to the ordering temper- 
ature Tjv despite its small size. As will be shown in a 
forthcoming publicationE2l the interchain coupling has a 
strong effect on the excitation intensity for all tempera- 
tures T < 50 K and there is considerably less magnetic 
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intensity in Sn(Q) as well as in S(Q) near Q c = 1 than 
in uncoupled chains. We may also imagine as T/v is ap- 
proached that spectral weight is transferred from the ID 
point to the 3D ordering wave-vector. 



B. The Dispersion Relation 

The spectrum of excitations of CSN1CI3 was measured 
for T = 12, 25 and 50 K over a wide range of wave- 
vector transfers using the MARI spectrometer with the 
incident beam perpendicular to the c-axis. The neutron 
scattering was measured for energy transfers of up to 75% 
of the incident energy. For a ID magnet, the scattering 
S(Q,uj) was obtained as a function of Q c and u>, where 
Q c is the component of the wave- vector transfer along 
the chain, while the in-plane wave-vector Q a is a known 
function of Q c and of the energy transfer. 
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FIG. 6. Neutron scattering spectra at T = 25 K for three 
different Q c . Open circles show the background as explained 
in the text. The solid line is the antisymmetrized Lorentzian 
weighted with the Bose factor. 



For measurements below T < 25 K, where the excita- 
tions may exhibit a dispersion perpendicular to the chain 
direction, only data only from the central detector banks 
were taken to limit the out-of-plane wave-vectors sam- 
pled in the experiment. Figs. ^ and [7] show the observed 
neutron scattering intensity for three wave- vector trans- 
fers along the chain at T — 25 and 50 K. With increasing 



temperature, the intensity for Q c = 1 decreases, moves 
to higher energies and broadens, as discussed in the pre- 
vious section. 

Our previous study has shown that,the excitation at 
Qc = 1 remained resonant up to 70 K.c3 Figs. [] and fj] fur- 
ther show that the excitations remain resonant through- 
out the zone and that they are weakly dispersive even at 
50 K. The scattering was again htted to the Lorentzian 
of Eq. U but now taking the line-shape of the quasi-elastic 
peak for the energy resolution. Least-square fits gave the 
parameters of Eq. A, e(Q) and T (see in Fig. ||). Data 
with Q c > 1.3 was not analyzed because it was contami- 
nated with phonon scattering. 
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FIG. 7. Neutron scattering spectra at T = 50 K for three 
different Q c . Open circles show the background as explained 
in the text. The solid line is the antisymmetrized Lorentzian 
weighted with the Bose factor. 



The non-magnetic background was estimated from the 
scattering observed in detectors at high scattering angles 
where the wave- vector transfer |Q| is large, the magnetic 
form factor is low and the spectrum is dominated by the 
phonon scattering which is proportional to |Q| 2 . The 
non-magnetic background for each scan was obtained by 
scaling the intensity according to |Q| 2 and accounting for 
background scattering that is independent of |Q|. The 
latter was estimated from the energy gain side of the 
spectrum at low temperatures. The background is shown 
in Figs. H and ^|by open circles. Assuming that the whole 
intensity at high scattering angle scales only with |Q| 2 , 
or that part of the scattering scales with |Q| 4 leads to 
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even lower background estimates for the low-angle data. 
The results show that the background is weakly depen- 
dent on energy transfer and therefore may be accurately 
subtracted. 

The dispersion of the excitations decreases with in- 
creasing temperature as shown in Fig. ||a. While the ex- 
citation energy increases with increasing temperature for 
antiferromagnetic fluctuations, it decreases for Q c = 0.25 
and 0.5. Due to the lack of a theoretical model for the 
dynamic structure factor for the whole magnetic zone the 
dispersion was phcnomcnologically described by 

u 2 (Qc,T) = 

A 2 (T) + v 2 (T) sin 2 (7rQ c ) + a 2 (T) cos 2 (^) . (8) 

Here A(T) is the Haldane gap at Q c = 1, v s (T) is the 
spin velocity determining the increase of the excitation 
energy as \Q C — 1| increases, and a(T) takes account of 
the two-particle peak at low T as Q c — > (a 2 = 3), but 
in practice a is used as a free parameter to account for 
the asymmetry about Q c = 0.5. For T < 25 K 3D ef- 
fects were taken into account using a RPA approximation 
which introduces a small dispersion perpendicular to the 
chain axis.til The dispersion then reads 

oj 2 (Q) = lu(Q c ) ■ [u(Q e ) + 2 • J'(Q a , Q b , 0) ■ 5(Q)] . (9) 

where J'(Q a , Qb, 0) is the Fourier transform of the inter- 
chain coupling and S(Q) = f S(Q,<j)du) is the structure 
factor. The dispersion perpendicular to the chain is most 
pronounced at Q c = 1 where 5(Q) is maximum and is 
small at Q c = 0.5. 
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FIG. 8. (a) Excitation energy e(Q c ) as a function of 
wave-vector and temperature. The lines are fits described 
in the text, (b) Lorentzian half width I\ 



The Haldane gap or resonance energy A, the spin ve- 
locity v s and a were obtained for the different tempera- 
tures by fitting Eq. ^ to the Q-dependence of the excita- 
tion energies. The results for A(T) are shown in F-igs. ||a 
and H and are consistent with previous results. E2l The 
results for v s (T) and a 2 (T) are shown in Fig. The 
spin velocity v s is approximately constant below 15 K 
at 5.70(7) meV and above 15 K it decreases with tem- 
perature reaching 3.4(1) meV at 50 K. The asymmetry 
parameter a 2 is 6 meV 2 at T = 6 K, decreases rapidly 
with increasing temperature and becomes negative above 
about 12 K. Since it is merely a phenomcnological pa- 
rameter, this is not unphysical but reflects the reversal of 
the asymmetry about Q c — 0.5 as the high temperature 
limit is approached. 

The excitations strongly broaden with increasing tem- 
perature (Fig. ||b). Below T — 25 K, the width is ap- 
proximately independent of wave- vector. For Q c > 1.2 
the width begins to increase and we have evidence that 
the scattering there is contaminated with phonons and 
so these data were not further analyzed. For T — 50 K 
the excitation width increases from Q c = 0.2 to Q c = 1 - 
an indication of a change in its Q c -dependence with in- 
creasing temperature. We shall see that this is expected 
at high temperatures. The widths shown in Fig. ||b are 
slightly smaller relative to those in Fig. [|, presumably 
because the resolution width was overestimated in our 
analysis. 




10 20 30 40 50 

T(K) 

FIG. 9. The upper panel shows the spin velocity v s and the 
lower panel a 2 as a function of temperature. The solid lines 
are guides to the eye. 



C. Structure factor and correlation length 

The structure factor S(Q) was determined numerically 
from the measured neutron spectra at positive energy 
transfer after subtraction of the flat background. The 
intensity at negative energy transfers was estimated using 
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detailed balance. The structure factor was corrected for 
the wave- vector dependence of the magnetic form factor 
and for the effect of the 3D interactions by using Eq. || 
to derive the structure factor for the fluctuations at the 
ID point: 



^/vl +a 2 /4 



S(Q ) 



u(Q c ) 



S(Q). 



(10) 



This can be thought of as the structure factor of a single 
chain embedded in the system. 



O 

(XI 
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FIG. 10. Structure factor S(Q C ) along (0.81, 0.81, Q c ) for 
T = 12, 25 and 50 K after an allowance for 3D effects (Eq. 0). 
S(Q C ) was corrected for background and magnetic form fac- 
tor. The solid lines correspond to a fit to a square-root 
Lorentzian (Eq. |l^ ) as explained in the text. 



The effective structure factor for the ID chain (Fig. |10| ) 
shows strong antiferromagnetic fluctuations at Q c = 1 
that persist to T ~ 2 J. They decline with tempera- 
ture while S(Q C ) at Q c — 0.25 and 0.5 increases. For a 
ID magnet with a magnetic excitation spectrum that is 
dominated by a single mode, the Q c -dependence of the 
structure factor S aa (Q c ), a = x, y, z, can be expressed 

aska 



S aa (Q c ) = - 



4 < H > (1 - costtQ c ) 
3 L Tiu> aa (Q c ) 



(11) 



Here <TL >/L is the ground state energy per spin and 
uj aa (Q c ) is the single mode energy. This approximation 
is called the single mode approximation (SMA). Tak- 
ing the energy from Eq. ||, the structure factor close to 
the antiferromagnetic point behaves like a square-root 
Lorentzian: 



S(Q C ) oc 



where we expect 



(12) 



(13) 



and where £ is the correlation length. £ was determined 
from S(Q C ) for 0.7 < Q c < 1.3 by fitting Eq. || to the 
data shown in Fig. |l0|. 

The temperature dependence of the correlation length 
£ is shown in Fig. [ll] on a semi-logarithmic plot. At 50 K 
it is only slightly longer than 1 site which means that 
the spatial correlations are extremely short ranged. It 
is surprising to have resonant and dispersive modes in a 
spin system with such a short correlation length. With 
decreasing temperature, the correlation length increases 
to about 4 sites below 8 K, as previously reported.EII 
These results can be compared to the correlation length 
obtained by Kim et a/.ES from numerical quantum Monte 
Carlo calculations in which effective spin S chains were 
realized by ferromagnetically coupling.-a = 2S antifer- 
romagnetic spin chains with S = 1/2. cB Our measured 
correlation length for temperatures above 12 K is in ex- 
cellent agreement with the numerical results, and only 
below 12 K where 3D correlations are most important is 
the correlation length shorter than the numerical predic- 
tion for uncoupled spin-1 chains. That it becomes shorter 
is associated with the fact that we are at the non-critical 
ID wave- vector. 
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FIG. 11. Correlation length £(T) deduced from S(Q C ) and 
shown as a function of the reduced temperature T / ( JS(S+1)) 
on a semi-logarithmic plot. The solid line is the prediction of 
the correlat*ein length obtained from a quantum Monte Carlo 
calculation.!^! 



For a ID magnet in which the magnetic response is 
dominated by single mode excitations, the SMA (Eq. ||) 
predicts that £ is given by Eq. [l3|. In this experiment, 
all four quantities were independently determined and 
the relationship can be tested experimentally. Fig. [l2| 
compares the correlation length £ estimated from the in- 
dependence of the structure factor and the correlation 
length £ ca i c given by Eq. O. The agreement is excellent 
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for all temperatures. This demonstrates that SMA is 
valid up to T — 50 K and suggests that the observed 
excitations can be regarded as short-lived single particle 
excitations. 
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^ cal =(v s 2 + a 2 /4) 1,2 /A 

FIG. 12. The horizontal axis represents the correlation 
length £ ca i calculated from the experimentally determined A, 
v a and a for temperatures between T — 6 and 50 K. The ver- 
tical axis shows the correlation length £ determined from the 
wave- vector dependence of the structure factor S(Q C ). The 
solid line is £ ca i = £. 

The Hohenberg-Brinkman sum-ruleEl predicts that the 
first energy moment F(Q) = (1uiujS(Q,uj) of a ID 
magnet with isotropic nearest-neighbor exchange J is 
given by 

F(Q c ) = -^^^(l-cos(7rQ c )). (14) 
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Here Q c is the wave- vector transfer along the spin chain 
and < H >= 3J JT < S" S° +1 > is the expectation value 
of the Hamiltonian. 

The first energy moment F{Q C ) was determined nu- 
merically from the measured neutron scattering spectrum 
after subtraction of the background and accounting for 
the magnetic form factor, and is shown in Fig. O for 
three different temperatures. The dependence on Q c of 
the first moment, F(Q C ), is in excellent agreement with 
the sum-rule for all momenta and temperatures. The fit 
yields the internal energy < TL > in arbitrary units and 
its temperature dependence- and it will be discussed in 
a forthcoming publicationEj The data for < TL > are 
shown in Table 

IV. HIGH-TEMPERATURE PARAMAGNETIC 
SCATTERING 

At infinite temperature, there are no correlations be- 
tween the spins and the instantaneous structure factor 
S(Q) is a constant. The first moment of the spectrum 
vanishes for all Q because the Bose factor weighs the ex- 
citations for neutron energy gain and loss equally. The 
second moment, however, is predicted to vanish at the 
ferromagnetic momenta and reach a maximum at anti- 
ferromagnetic momenta, irrespective of the spin quantum 
number or whether the exchange is antiferromagnetic or 
ferromagnetic. E3E£| Chain systems with nearest-neighbor 
interactions are ideal systems to test this behavior, be- 
cause the wave- vector dependence of the second moment 
has a simple form: 

<lu 2 >cx J 2 [l-cos(7rQ c )], (15) 

where Q c is the wave- vector transfer along the chain. 

The scattering was measured at T = 200 K ~ 10J for 
a wide range of wave-vector transfers using the MARI 
spectrometer. Fig. |l4| shows the scattering surface pro- 
jected onto the (Q c , lj) plane for neutron energy loss. Be- 
cause the effect of the interchain interaction is negligible 
at T = 200 K, the data from all the detectors were used 
to improve the statistics. Fig. [ll| shows that the scatter- 
ing consists of broad scattering and that there is no sign 
of the resonant excitations seen at lower temperatures. 
The broad dome of spectral weight is centered at Q c = 1 
for the high energy fluctuations whereas it was located 
at Q c = 0.5 at low temperatures. 

Fig. [l5| shows constant-Qc scans for different wave- 
vector. The upper bound of the scattering has a maxi- 
mum of 10 meV at Q c = 1 and decreases with increasing 

|Qc-l|. 



FIG. 13. First energy moment F{Q C )= f°° dui uj S(Q, lj) for 
T = 12, 25 and 50 K. The solid line is a°°fit of the Hohen- 



berg-Brinkman sum rule (Eq. l 
to the internal energy < TL >J 



to a parameter proportional 
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FIG. 14. The intensity measured at T = 200 K using 
MARL The data were measured with an incident energy 
Ei = 20 meV and the intensity in all three detector banks was 
added and smoothed. The intensity is indicated by the bar on 
the right-hand side. The black trapezoidal areas are an arti- 
fact and a result of a MARI detector gap in the 26 scattering 
angle and the data smoothing algorithm. The quasi-elastic 
scattering below ~ 1 meV should be ignored. 



The scattering intensity observed at the antiferromag- 
netic point Q c — 1 is not consistent with the predifc 
tion of a coupled-mode theory at infinite temperature,E3 
shown as a dashed line in Fig. 15. Instead the scatter- 
ing spectrum rias the form of a Gaussian as predicted 
by de GennesH (solid line). The difference between the 



coupled-mode theory and the experiment may be because 
the experiment was performed at a finite temperature but 
more likely because the theory is unsatisfactory. 

The structure factor S(Q C ) and the first energy mo- 



16| in the same units as for 
and Il0|). They were cal- 



ment F(Q C ) are shown in Fig 
the lower temperatures (Figs 
culated numerically from the observed neutron scattering 
at positive neutron energy transfers after subtraction of 
a flat background that matched the scattering at high en- 
ergy transfers. The intensity at negative energy transfers 
was estimated using detailed balance. S(Q C ) and F(Q C ) 
were corrected for the magnetic form factor. S(Q C = 1) 
is now only about twice S(Q C = 0.25), showing that there 
are weak spin correlations, but these are considerably re- 
duced compared to lower temperatures. The solid line 
is the prediction of a high-temperature series expansion 
Eq. 28 ) which will be discussed in more detail in Section 




0.5 
0.0 

5 

Energy (meV) 

FIG. 15. Neutron spectra at T = 200 K for four different 
Q c . The peak centered at zero energy transfer arises from 
quasi-elastic incoherent scattering from the sample and the 
sample holder. The open circles show the background as ex- 
plained in the text. The dashed and solid lines are the psedic- 
tion of a coupled-mode theory at infinite temperaturej£j and 
the de-Gennes Gaussian fitted to the data, respectively. The 
instrumental resolution (FWHM) is shown by a bar. 



The Qc-dependence of the first moment F(Q C ) is 
well described with the sum rule of Hohenberg and 
BrinkmarS, confirming that the scattering is magnetic 
and that contamination with phonon scattering is not 
significant. The maximum of F(Q C ) is ~ 6.5 times lower 
at 200 K than the maximum at 12 K, showing that al- 
though the first moment is not zero, the system is close 
to the high-temperature limit. 

It is clear from Figs. |IJ and [l5| that the second moment 
has a maximum at the antiferromagnetic point Q c = 1 as 
predicted by the theory of de Gennes. It will be shown 
in the section |y| that there is quantitative agreement 
between experiment and theory if the scattering is ex- 
panded to first order in 1/ksT. 



10 



=3 

o 

o 
en 



=3 

3, 6 
^ 2 



a) 
















CsNiCI 3 


T = 200 K 


1 1 1 1 

b) 


1 1 1 — 








T = 200 K 



0.2 0.3 0.4 0.5 0.6 0.7 

Q c (r.l.u.) 



0.9 1.0 



FIG. 16. (a) S(Q C ) as a function of wave-vector transfer 
Q c at T — 200 K in the same units as in Fig. [l^. The solid 
line is the prediction of the high-temperature series expan- 
sion (Eq. normalized to the data, (b) Observed first 
energy moment F(Q C ). The solid line is the fit of the Hohen- 
berg-Brinkman sum rule, ~f <r v > (1 — cos(7rQ c )). 



V. DISCUSSION 

A. Life-time of Haldane excitation for T < A 

The excitations of an antiferromagnetic isolated spin-1 
chain in its quantum-disordered phase are well-defined, 
gapped spin-1 excitations (Haldane excitations) which 
satisfy a relativistic dispersion relation, e(q) = (A 2 + 
\Q — 1|. At zer^temperature where 



u 2 g 2 ) 1/2 , where q 



the spin correlation length is 6 sites,E£lthe neutron injects 
an excitation in the Haldane band which travels with a 
group velocity v 2 q/e(q), and its life-time is infinitely long. 
At a finite temperature T < A = 0.41J, thermal activa- 
tion creates a dilute gas of Haldane excitations. They 
limit the life-time of the injected excitation through col- 
lisions, and this gives it a finite energy E3 

The scattering matrix for a two-particle collision has 
a superuniversal form in ID which does not depend on 
any microscopic parameter, and for T -C A, it predicts 
that a two-particle state ({qi, <j{\, {q2, 02}) transforms 
into a final state ({91, 02}, {(72, 0i}) 5 where a is the_spin 
quantum number, thus basically swapping the spin.u Al- 
ternatively in ID we may think that the particle created 
by the neutron, cti}, collides with one of the gap par- 
ticles {92)02} and, like billiard balls, the momentum of 
the first is transmitted to the second, to give {92, 01} 
and {(71,02}. This result holds for any_form of ID dis- 
persion relation. Damle and SachdevEfl calculated the 
temperature dependence of the line-shape and width by 
approximating the relativistic dispersion with a classi- 



cal dispersion, e(q) — A + and when T <C A, the 

dynamic structure factor has a simple scaling form 



(16) 



where 7 is the inverse of the time between collisions with 
other excitations and A is a scaling constant. The scaling 
function $(z) was calculated numerically and found to be 
close to a Lorentzian form: 



ira 



2(a 2 + z 2 ) 



(17) 



with a = 0.71. Thus the theory predicts that the ob- 
served half- width is T — a ■ 7. The inverse collision time 
is proportional to the density of the excitations times 
their root mean square velocity, namely 



7(T) = — — exp 



A 



(18) 



the same as for classical gas collisions. The experimen- 
tally observed excitation T width in CsNiCl3 is apprecia- 
bly lower than this prediction as shown in Fig. for all 
T<A~11K. 

We have modified the calculation of the life-time of 
the excitations by starting from the proper relativistic 
dispersion instead oLthe classical approximation used by 
Damle and SachdevO We have also applied a finite lat- 
tice momentum cut off. The density of the quasi-particle 
gas is then given by 



Prel(T) 



dq 
2^ 



exp 



V^ 2 



v 2 q 2 



knT 



(19) 



The mean square velocity of the quasi-particles is given 
by 



i(T) = 



Prcl(T) J-n 27T 



d( i f 2 n2 

— (v s q) exp 



V^ 2 



v 2 q 2 



knT 



(20) 



These two integrals are not analytically solvable but were 
calculated numerically to obtain 7 re i(T') = v rc \(T) p rc \(T) . 
7 TO i is larger than 7 because the excitation energy of rela- 
tivistic particles is lower at large momenta than for clas- 
sical particles with the same A and v s . At the same 
temperature more excitations are then populated so that 
collisions are more frequent. 

r rc i = a ■ 7 rc i is shown in Fig. ^ as dashed line, and 
in the expected range of applicability the agreement be- 
tween this calculation and the experiment is good for 
T > 8 K. This agreement is evidence for the validity of 
the model by Damle and Sachdev that describes an an- 
tiferromagnetic spin-_L chain at T < A as a dilute gas of 
Haldane excitations .EJ Our improved calculation merely 
helps the theory to remain valid to the relatively high 
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temperatures of our experiment. Below T < 8 K the 
observed excitation widths are higher than the calcula- 
tions, most probably due to the interchain spin corre- 
lations which soften the Haldane gap near the ordering 
wave-vector (0.33,0.33,1), so that there is an increased 
thermal population of excitations. 

B. Temperature dependence of solitonic excitations 
in quantum-disordered phase 

An antiferromagnetic spin-1 Heisenberg chain does not 
order at any finite temperature and the spin-spin cor- 
relations decay exponentially with distance. However, 
there is a form of hidden spin order given by a non-local 
string ordet, parameter first proposed by Den Nijs and 
Rommelse.E3 It is defined as 

0: tYing (H) = limi^i^ < of. > n , (21) 

where 

^ = -Sfe^(m ( 22 ) 

V i=i+i ) 

Kennedy and TasakiEl showed the hidden spin order can 
be regarded as consequences of the breaking of the Zi x Zi 
symmetry. It was shown that 0" tling (H) — 0.374 for 
an antiferromagnqtip|Ghain with nearest-neighbor Heisen- 
berg interactionsBacJ The hidden spin order can be 
pictured, in the frame work of the valence solid bond 
(VBS)Jl^l where the 5 — 1 spins are obtained by sym- 
metrizing two 5=1/2 spins on one site, and each of the 
5 = 1/2 spins forms a coherent singlet state with an- 
other 5 = 1/2 spin from a neighboring site. A possible 
VBS spin configuration for the ground state is shown in 
Fig. [17](a) together with its configuration in the spin-1 
picture. The S z = 1 and S z = — 1 states alternate with 
a random number of S z = states in between, and the 
ground state can be viewed as diluted spin-1 Neel or- 
der, which is the string order parameter defined above. 
In VBS (D° ling (H) = f . This value is higher than that 
for the Heisenberg model because quantum fluctuations 
reduce the long-range hidden order in the latter. 

A scattered neutron excites a spin- 1 /2 p air into a sym- 
metric triplet state, as shown in Fig. |l7j(b). The corre- 
sponding spin configuration in the spin-1 picture shows 
that the dilute spin-1 Neel order is now broken. Fath and 
Solyom have shown that such a triplet states are moving 
hidden domain walls and are coie-soliton excitations of 
the string-ordered ground state.u 

At finite temperature, the hidden non-local string or- 
der is destroyed because spin flips are allowed through 
thermal population of the Haldane band that change the 
instantaneous jnagnetization of the chain. Yamamoto 
and Miyashitac3 introduced a local (short-range) hidden 
order parameter 



SR N(+-) + N(-+)+N(++)+N{—) [ ' 

where V(±+) and iV(±±) are the total number of near- 
est spin pairs of (±+) and (++) type, respectively. The 
operator is applied to a spin wave-function, for example 
(.. — OH — + — OH — + 0— ..), where all spin projection 

are omitted, leaving (.. — | 1 1 1 — ..). The 

expectation value < Osr > is 1 in the absence of (++) 
and ( ) pairs which can regarded as domain walls. 
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FIG. 17. (a) VBS ground state where the spin-1 (large 
circles) can be pictured as two symmetrized spin- 1/2 (small 
circles within the large circles). Two spin- 1/2 from neigh- 
boring sites form singlet bonds (straight lines between two 
spin- 1/2). The picture shows one possible spin- 1/2 configura- 
tion in the ground state and the corresponding spin-1 config- 
uration showing the hidden spin string order, (b) VBS state 
after a neutron has transferred S z = — 1 to a spin-1, creating 
a symmetrized triplet state between two spin-1/2 on neigh- 
boring sites (the cross between two spin-1/2). This breaks 
the hidden spin-1 order. 

While the non-local string order collapses at any finite 
temperature, the local string-order decreases only slowly 
with increasing temperature,E3 and reduces only by a fac- 
tor of 2.7 from T = to 50 K. We note that S{Q C = 1) 
grows by 1.5 between 50 K and 12 K and is not too 
much less than the factor of 2 growth of the local string 
order .c3 In the light of the persistence of the hidden spin 
order to high temperatures it is tempting to interpret the 
observed excitations as moving domain walls of the spin 
string order, and it is less surprising that they are reso- 
nant at such high temperatures. More theoretical work 
is needed to clarify this point, and to describe the energy 
and the life-time of these high-temperature solitons. 

C. High-temperature paramagnetic scattering 

In the high-temperature limit J(3 = J /ksT <C 1, the 
second moment of the dynamic structure factor is given 
byM 
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< or >= 



-T.i.eMiQiRi- Rj)) < m,H],s<?] > 



(3 < S%S a _ Q > 



(24) 



The thermal average of the spin-spin correlations is 



< SfS* >= Tr (exp(-/3H)S?S?) /Tr (exp(-/?H)) 



To first order in /3, 



and 



< SfS* >=<5 Q , 5 -S(S + 1) 



(25) 



(26) 



<SfS d j >=-S a , s l3\-S(S + l)j JiRt-Rj). (27) 

Please note that we use a different definition of J that 
that used by Loveseyo. The Fourier transform of the 
spatial spin-spin correlations < SfS^ > is 



1 



< S%S! Q >« -S(S +1)1- -/3S(S + 1) J(Q) 



1 



(28) 



For a ID magnet with a nearest-neighbor exchange J, 
J(Q) — 2Jcos(7rQ c ) and the second moment becomes 



< t*r >= - 



S(S+1) 



3 1- yS(S+l)J(Q c ) 



J 2 (1-cos(ttQ c )) 
(29) 



This result can be directly compared with the mea- 
surements at T — 200 K. The calculated S(Q C ) given 
by Eq. |28| with J = 2.28 meV is in excellent agreement 
with the experiment at T — 200 K (Fig. |l(|) if an overall 
scaling factor is fitted to the experimental data. 
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FIG. 18. Square root of the second energy moment fl(Q c ) 
for several temperatures as a function of wave- vector transfer 
Qc. The solid line is the theoretical prediction given by Eq. 
The dashed line is the prediction for infinite temperature. 



TABLE I. For neutron energy loss, the integrated intensity Ro(Q) = dcoS(Q, cu) and the first and second moment 
< w > (Q) — J °° duiS(Q, uo)uj / Ro(Q) and < cu 2 > (Q) = J Q duiS(Q, w)uj 2 / Ro(Q) of the observed spectra. Ro is in arbitrary 
units (a.u.) and < uj > and < ui 2 > are in units of meV and meV 2 , respectively. Columns 2 to 4 are derived from MARI 
data with a varying wave-vector transfer perpendicular to the chain. The final columns are derived from the RITA data with 
a different normalization. 
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10 
















1404(32) 


1.95(0.06) 


2.36(0.05) 


12 


4.3(1.0) 


5.8(0.6) 


5.9(0.4) 


20.84(0.66) 


2.89(0.6) 


3.58(0.22) 


-26(0.5) 


1233(45) 


2.17(0.13) 


2.78(0.13) 


15 
















1055(42) 


2.17(0.16) 


2.52(0.23) 


20 
















939(39) 


3.1(0.2) 


3.62(0.13) 


25 


6.1(1.0) 


5.0(0.4) 


5.4(0.3) 


16.55(0.48) 


3.72(0.2) 


4.37(0.24) 


-21.3(0.2) 


813(37) 


3.04(0.22) 


3.66(0.15) 


30 
















729(36) 


3.0(0.23) 


3.48(0.18) 


35 
















689(35) 


3.37(0.25) 


4.03(0.14) 


50 


5.3(0.9) 


4.4(0.3) 


4.9(0.3) 


10.89(0.26) 


4.03(0.12) 


4.49(0.1) 


-14.2(0.2) 


641(29) 


3.75(0.23) 


4.43(0.09) 


70 
















571(28) 


3.78(0.26) 


4.43(0.1) 


200 


5.2(0.9) 


3.5(0.2) 


4.0(0.2) 


7.5(0.19) 


4.46(0.14) 


4.95(0.12) 


-3.8(0.1) 
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A more sensitive test of the theory is a quantita- 
tive comparison of the predicted and measured second 
energy moment. The square root of the second mo- 
ment Q =< lu 2 > x / 2 was determined numerically from 
the measured spectra after the usual corrections were 
made. Fig. |l8| shows that £1{Q C ) has a maximum which 
approaches Q c = 1 at high temperatures. The solid 
line shows the predicted fl{Q c ) according to Eq. ^9] at 
T = 200 K and using J = 2.28 meV. There is good 
quantitative agreement between theory and experiment 
considering that the solid line is not a fit but a prediction 
based only on the temperature and the exchange constant 
along the chain. Q(Q C ) for infinite temperatures (3 = 
is shown by the dashed curve for comparison. 

Fig. [l^ also shows that the maximum in fl(Q c ) changes 
from Q c — 1 to Q c — 0.5 on cooling. At T = 12 K, the 
maximum is where the dispersion has its largest energy. 
The moments of the spin spectrum are summarized in 
Table |. Here the integrals are done only over the positive 
energy response. 

VI. CONCLUSIONS 

In summary, the evolution of the dynamic spin struc- 
ture factor S{Q c ,u) of CsNiCkj has been surveyed be- 
tween T = 5 and 200 K over a wide range of wave- vector 
transfers Q c along the chain. 

The temperature dependence of the correlation length 
£(T) was determined from the Q c -dependence of the 
equal-time structure factor and the results are in agree- 
ment with quantum Monte Carlo calculations except for 
T < 10 K, close to the ordering temperature Tm- At 
these temperatures and at the non-critical ID wave- 
vector (0.81, 0.81, 1) the correlation length is shorter, sur- 
prisingly, than is predicted for isolated chains, and the 
gap is higher than expected. The intensities of the ex- 
citations do not scale with the 1 / w-law which is typical 
for antiferromagnets. These results are in sharp contrast 
with the predictions of the random phase approximation 
and with the single mode approximation. 

The magnetic excitations remain resonant and disper- 
sive up to T — 50 K ~ 2J, a temperature comparable 
with the spin-band width. They remain resonant because 
the energy of the spin excitations renormalizes upward 
more rapidly than their damping. The temperature de- 
pendence of the dispersion along the chain direction was 
measured, and the excitation energy A(T) and the spin 
velocity v s (T) were determined for 6 K < T < 50 K. 

,/v 2 +a 2 /4 

We found that the relation £ = v " A , which is pre- 
dicted by the single mode approximation, remains valid 
in the entire temperature region. This suggests that the 
excitation spectrum around the antiferromagnetic point 
consists of single-particle excitations even at T — 50 K 
where the excitations are broad. It also hints at the exis- 
tence of a hidden spin order which carries the propagating 
excitations as predicted by numerical calculations. The 



observed excitations may be regarded as domain walls 
(solitons) of the hidden string order, whose energy scale 
is J rather than A. 

At T = 200 K the scattering is well described as 
paramagnetic scattering and extends to ~ 10 meV en- 
ergy transfer at the antiferromagnetic point. The upper 
boundary of the paramagnetic scattering exhibits a broad 
maximum at the antiferromagnetic point. 

The structure factor and the second energy moment 
are consistent with the first-order high-temperature se- 
ries expansion. These results call for an extension to 
higher temperatures of the theory for gapped weakly cou- 
pled spin chains. 

ACKNOWLEDGMENTS 

We would like to thank Dr. O. Petrenko for his assis- 
tance with experiments at ISIS and Dr. P. Santini and 
Dr. D. A. Tennant for enlightening discussions. We are 
grateful to K. N. Clausen and the staff of Ris0 National 
Laboratory for their expert assistance and for the use of 
the excellent neutron facilities that made this research 
possible, regrettably no longer available with the shut 
down of DRB. Financial support for the experiments was 
provided by the EPSRC, by the EU through its Large In- 
stallations Program and by the British Council-National 
Research Council Canada Program. ORNL is managed 
for the U.S. D.O.E. by UT-Battclle, LLC, under contract 
no. DE-AC05-00OR22725. One of the authors (M. K.) 
was supported by the Swiss National Science Foundation 
under Contract No. 83EU-053223. 



1 L. D. Faddeev and L. A. Takhatajan, Phys. Lett. A 85, 
375 (1981). 

2 F. D. M. Haldane, Phys. Rev. Lett. 50, 1153 (1983). 

3 F. D. M. Haldane and M. R. Zirnbauer, Phys. Rev. Lett. 
71, 4055 (1993). 

4 R. Botet and R. Julien, Phys. Rev. B 27, 613 (1983). 

5 M. P. Nightingale and H. W. J. Blote, Phys. Rev. B 33, 
659 (1986). 

6 W. J. L. Buyers, R. M. Morra, R. L. Armstrong, M. J. 
Hogan, P. Gerlach, and K. Hirakawa, Phys. Rev. Lett. 56, 
371 (1986). 

7 M. Steiner, K.Kakurai, J. K. Kjems, D. Petitgrand, and R. 
Pynn, J. Appl. Phys. 61, 3953 (1987). 

8 J. P. Renard, L. P. Regnault, and M. Verdaguer, Journal 
de Physique 49, 1425 (1988). 

9 G. Fath and J. Solyom, J. Phys.: Condens. Matter 5, 8983 
(1993). 

10 D. A. Tennant, T. G. Perring, R. A. Cowley, and S. E. 
Nagler, Phys. Rev. Lett. 70, 4003 (1993). 

11 W. J. L. Buyers, Z. Tun, A. Harrison, J. A. Rayne, and 
R. M. Nicklow, Physica 180&181B, 222 (1992). 



14 



12 I. A. Zaliznyak, L.-P. Regnault, and D. Petitgrand, Phys. 
Rev. B 50, 15824 (1994). 

13 H. Mutka, C. Payen, P. Molinie, J. L. Soubeyroux, P. 
Colombet, and A. D. Taylor, Phys. Rev. Lett. 67, 497 
(1991). 

14 J. Darriet and L.-P. Regnault, Solid State Commun. 86, 
409 (1993). 

15 H. Mutka, C. Payen, P. Molinie, and R. S. Eccleston, Phys- 
ica B 213&214, 170 (1995). 

16 T. Sakaguchi, K. Kakurai, T. Yokoo, and J. Akimitsu, J. 
Phys. Soc. Jpn. 65, 3025 (1996). 

17 J. P. Renard, M. Verdaguer, L.-P. Regnault, W. A. C. Erke- 
lens, J. Rossat-Mignod, J. Ribas, W. G. Stirling, and C. 
Vettier, J. Appl. Phys. 63, 3538 (1988). 

18 S. Ma, D. H. Reich, C. Broholm, B. J. Sternlieb, and R. W. 
Erwin, Phys. Rev. B 51, 3289 (1995). 

19 R. M. Morra, W. J. L. Buyers, R. L. Armstrong, and K. 
Hirakawa, Phys. Rev. B 38, 543 (1988). 

20 M. Kenzelmann, R. A. Cowley, W. J. L. Buyers, and D. F. 
McMorrow, Phys. Rev. B 63, 134417 (2001). 

21 A. A. Katori, Y. Ajiro, T. Asano, and T. Goto, J. Phys. 
Soc. Jpn. 64, 3038 (1995). 

22 N. Achiwa, J. Phys. Soc. Jpn. 27, 561 (1969). 

23 K. Lefmann, D. F. McMorrow, H. M. Ronnow, K. Nielsen, 
K. N. Clausen, B. Lake, and G. Aeppli, Physica 283B, 343 
(2000). 

24 K. Damle and S. Sachdev, Phys. Rev. B 57, 8307 (1998). 

25 M. J. Cooper and R. Nathans, Acta Crys. 23, 357 (1967). 

26 T. Jolicceur and O. Golinelli, Phys. Rev. B 50, 9265 (1994). 

27 M. Kenzelmann, R. A. Cow ley, W. J. L. Buye rs, Z. Tun, 
R. Coldea, and M. Enderle, |sond-mat/01 12152 . 

28 M. Kenzelmann, R. A. Cowley, W. J. L. Buyers, R. Coldea, 
J. S. Gardner, M. Enderle, D. F. McMorrow, and S. M. 
Bennington, Phys. Rev. Lett. 87, 017201 (2001). 

29 Y. J. Kim, M. Greven, U. J. Wiese, and R. J. Birgeneau, 
Eur. Phys. J. B 4, 291 (1998). 

30 M. Kenzelmann and P. Santini, Temperature dependence 
of single particle excitations in a 5=1 chain: exact di- 
agonalization calculations compared to neutron scattering 
experiments, submitted to Phys. Rev. B. 

31 Z. Tun, W. J. L. Buyers, R. L. Armstrong, K. Hirakawa, 
and B. Briat, Phys. Rev. B 42, 4677 (1990). 

32 P. C. Hohenberg and W. F. Brinkman, Phys. Rev. B 10, 
128 (1974). 

33 P. G. Gennes, in Magnetism, Vol. 3 of G. T. Rado and H. 
Suhl, p. 115, (Academic, New York, 1960). 

34 S. W. Lovesey, Theory of neutron scattering from condensed 
matter (Clarendon Press, Oxford, 1984). 

35 S. W. Lovesey and E. Balcar, J. Phys.: Condens. Matter 
6, 1253 (1994). 

36 M. den Nijs and K. Rommelse, Phys. Rev. B 40, 4709 
(1989). 

37 T. Kennedy and H. Tasaki, Phys. Rev. B 45, 304 (1991). 

38 S. M. Girvin and D. P. Arovas, Phys. Scr. T 27, 156 (1989). 

39 Y. Hatsugai and M. Kohmoto, Phys. Rev. B 44, 11789 
(1991). 

40 S. R. White and D. A. Huse, Phys. Rev. B 48, 3844 (1993). 

41 I. Affleck, J. Phys.: Condens. Matter 1, 3047 (1989). 

42 S. Yamamoto and S. Miyashita, Phys. Rev. B 48, 9528 
(1993). 



15 



